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An impurity coupling to both spin and charge degrees of freedom is added to a periodic t-J 
chain such that its interaction with the bulk can be varied continuously without losing integrability. 
Ground state properties, impurity contributions to the susceptibilities and low temperature specific 
heat are studied as well as transport properties. The impurity phase-shifts are calculated to establish 
the existence of an impurity bound state in the holon sector. 



Oh! 



> 
(N 

NO 

(N 
On 
O 
NO 

On 
c3 



C 

O 
O 



Quantum fluctuations are known to play an important 
role in the physics of low dimensional strongly correlated 
electron systems: the low temperature properties of such 
systems in one spatial dimension have to be described in 
terms of a Luttinger liquid rather than a Fermi liquid. 
From an experimental point of view the transport prop- 
erties of these systems in the presence of boundaries and 
potential scatterers are of particular interest. Recently 
several attempts have been made to describe such a situa- 
tion: Using rcnormalization group techniques the trans- 
port properties of a ID interacting electron gas in the 
presence of a potential barrier have been studied by Kane 
and Fisher jy. Their surprising findings triggered fur- 
ther work using different techniques like boundary con- 
formal field theory Q and an exact solution by means 
of a mapping to the boundary sine-Gordon model 
In particular the low temperature properties of magnetic 
(Kondo) impurities in a Luttinger liquid ||,[| have been 
investigated in great detail. In the present work we will 
investigate the effects of a particular type of potential im- 
purity in a Luttinger liquid (where both spin-and charge 
degrees of freedom are gapless) by means of an exact so- 
lution through the Quantum Inverse Scattering Method 
(QISM) §. 

Attempts to study effects due to the presence of impu- 
rities in many-body quantum system in the framework of 
integrable models have a long successful history [^[dl2[. 
As far as lattice models are concerned the underlying 
principle in these exact solutions is the fact that the 
QISM allows for the introduction of certain "inhomo- 
geneities" into vertex models without spoiling integra- 
bility. The local vertices — so called ^-operators — are ob- 
jects depending on a complex valued spectral parameter 
acting on an auxiliary matrix space in addition to the 
quantum space of the model. They are solutions of a 
Yang-Baxter equation with an 1Z matrix which itself acts 
on two copies of the matrix space and depends on the 
difference of the corresponding spectral parameters only. 



This allows to to build families of vertex models with 
site-dependent shifts of the spectral parameters and even 
different quantum spaces on different sites. The first fact 
has been widely used in solving models for particles with 
an internal degree of freedom by means of the nested 
Bethe Ansatz fljl. The second approach has been first 
applied by Andrei and Johannesson to study the proper- 
ties of an S = \ Heisenberg chain with an additional site 
carrying spin S || (see also ||). 

In this letter we study the properties of the supersym- 
metric t-J model with one vertex replaced by an C oper- 
ator acting on a four-dimensional quantum space. This 
preserves the gl(2\l) supersymmetry of the model but 
at the same time lifts the restriction of no double occu- 
pancy present in the t-J model at the impurity site. The 
existence of a free parameter in the four-dimensional rep- 
resentations of the superalgebra |L4| allows to tune the 
coupling of the impurity to the host chain. As will be 
shown below, the present model allows to study some 
aspects of a more general situation than the ones men- 
tioned above: the impurity introduced here couples to 
both spin- and charge degrees of freedom of the bulk 
Luttinger liquid. The extension of our calculation to the 
case of many impurities is straightforward. 

The solution of the model is completely analogous to 
that of the pure t-J model |l5j : The transfer matrix gen- 
erating the hamiltonian and the other conserved quanti- 
ties is the trace of a product of the local C operators 
chosen as Cu = (A + ill) /(A + i) for the regular sites 
(LT is a graded permutation operator acting on the ten- 
sor product of the auxiliary and the quantum space) and 
£34 oc A — i(j + 1) + i£ on the site associated to the 
impurity. Written as a matrix in the three dimensional 
auxiliary space L reads 
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Here Q a = \fa- 



o*JuX 2 <t2 with a = ±, where 



+ (— ) corresponds to f (|) and where the Hubbard 
projection operators are given by X ab = \a)(b\ with 
a,b =|,J,,2,0. The hamiltonian is then given by the 
logarithmic derivative of the transfer matrix at spectral 
parameter A = 0. In general the form of this hamil- 
tonian is rather complicated and will be given elsewhere 
p6[ . In any case the precise form of the lattice (impurity) 
interactions is not essential as far as low-energy proper- 
ties are concerned: in the continuum limit only a small 
number of terms with scaling dimensions smaller than 
two will survive (taking the continuum limit and identi- 
fying the scaling dimensions of the composite operators 
at the impurity site is rather nontrivial though). Physi- 
cally the model describes an impurity with four allowed 
states (spin-up/down electrons, empty/doubly occupied 
site) that couples to two neighbouring t-J sites and also 
modifies the interaction between the t-J sites (see Fig. [l]). 
In the limiting cases a — > 0(oo) the hamiltonian simpli- 
fies essentially: In the first case the impurity acts as a 
ordinary t-J site in the ground state below half filling, for 
a — > oo the impurity site is doubly occupied and induces 
a twist in the boundary conditions of the host chain. 

The eigenstates of the hamiltonian for iVf (JVj_) elec- 
trons with spin f Q) are constructed by means of the 
nested algebraic Bethe Ansatz (NABA) leading to a sys- 
tem of algebraic equations for the spectral parameters Xj 
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The corresponding eigenvalues of the hamiltonian in the 
grand canonical ensemble are E = —fiN e — (H/2)(N-t — 

^;) + E?=ii/(Aj + i). 

The configuration of spectral parameters leading to the 
lowest energy state for given chemical potential and mag- 
netic field are found in complete analogy to the pure t-J 
chain: the ground state for finite H is described by two 
filled Fermi seas of A-A^- "strings" A± = A^ ± § with 
real A^ 1 -* associated with the holon excitations and real 
solutions Xj describing spin degrees of freedom. In the 
thermodynamic limit dressed energies e c and e s can be 
associated with the excitations of these objects. They 
are given in terms of coupled integral equations which 
are identical to those found for the chain without impu- 
rities Jl7| . In the resulting ground state energy the impu- 
rity contribution can be identified from its L-dependence 
which allows to compute the occupation, magnetization 
and susceptibilities of the impurity site. Analytical re- 
sults for these expectation values are available only in 
limiting cases close to half filling and for densities near 



or below the critical density n c related to the magnetic 
field by H = 4 sin 2 (7rn c /2) where the ground state is fer- 
romagnetic and only real Xj are present jl6| . For general 
values of band filling and magnetic field the magneti- 
zation and particle number on the impurity site can be 
determined by numerically solving a set of two coupled 
integral equations (see Fig. ^|): For a — * the impurity 
mimics the bulk behaviour as discussed above. For large 
a the impurity occupation number is close to 1 (2) for 
n e < n c (> n c ) leading to an enhanced (reduced) magne- 
tization at small (large) electron densities. The magnetic 
susceptibility of the impurity near half filling is found to 
be ~ Xbuik/ct- We note that the magnetization curves for 
sufficiently large a intersect the ones for small a. 

In addition to the ground state properties the Bethe 
Ansatz allows to study the finite temperature behaviour 
of the system. The thermodynamic Bethe Ansatz equa- 
tions for the system with impurity are the same as in the 
pure case [p"7f . Again the impurity contribution can be 
isolated in the free energy which can evaluated explicitly 
in the limit T — > oo and for H T using Takahashi's 
method uM . In the high temperature limit we find 
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giving the correct entropy in this limit. Note that the 
parameter a enters the leading term in this expansion in 
a trivial way only. 

For low temperatures T <C H we can determine the 
phase diagram of the system. Most interesting is the be- 
haviour at half filling where the impurity contribution to 
the specific heat is found to show a different temperature 
dependence than the one from the bulk: for 2/i > H > 4 
the system is ferromagnetic and the low T free energy is 
given by (we suppress the contribution from the ground 
state energy) 
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For smaller magnetic fields the thermodynamic equilib- 
rium state is not ferromagnetically ordered, the bulk free 
energy is Fbuik ~ —kLT 2 /(6v s ) with the spinon veloc- 
ity v s and the impurity contribution is Fi mp w — — T 3 / 2 
up to factors that cannot be caculated in closed form 

in general. Near H = H c this factor becomes ~ 

H)-i exp (i(ff/2 — jx + &(4 - Hfl^)). 

The effect of the impurities on the transport proper- 
ties of the system can be studied by calculating the spin- 
and charge stiffnesses from the finite size corrections to 
the ground state energy of the model subject to twisted 
boundary conditions |n|. Following the analysis in |2{j] 
we introduce twist angles 4> c and <p s affecting charge and 
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spin degrees of freedom, respectively. The leading term in 
the resulting shift of the ground state energy can then be 
written as AE(<fi c , 4> s ) = L~ 1 4> a D a p(u)4>i3 where charge 
(spin) stiffness are defined as flW = (L/2)djAE(4>,0) = 
D cc (D^ = (L/2)dlAE((f>,-24>)). The analysis of the 
finite-size corrections to the ground state energy yields 
the result that for the case of a single impurity the stiff- 
nesses are not modified to leading order in L . Hence, 
in spite of the presence of the impurity we find an in- 
finite dc-conductivity. This is completely different from 
the situation in the "weak-link" type potential impurity 
discussed in such a weak link drives the Luttinger 

liquid to a strong coupling fixed point characterized by a 
vanishing conductivity. We believe that the behaviour of 
the system considered here is related to its integrability 
and the absence of backscattering at the impurity. 

The transport properties of the system do change if 
one considers a finite density rii of impurity sites. In this 
situation the band filling can take values larger than 1 as 
the impurity sites allow for double occupancies. Compar- 
ing the charge stiffness to that of the pure t-J case one 
observes a reduction for densities just above the critical 
one. For larger band fillings the presence of the impuri- 
ties leads to an enhancement of the stiffness (see Fig. ||). 
This is easily understood: vanishes at half filling in 
the t-J chain. The impurities do allow double occupan- 
cies thereby enlarging the phase space for the electrons 
which leads to an increase in the stiffness. At large fill- 
ings the stiffness increases as a function of a as the aver- 
age occupation number of the impurity sites are close to 
doubly occupied which makes the movement of electrons 
between the t-J sites easier. In particular, the "absorp- 
tion" of particles by the impurity sites for a = oo leads to 
plateaus in the stiffness for n e < rii (where it vanishes) 
and for rii + (1 — rii)n c < n e < 2rii + (1 — m)n c . For 
other fillings the stiffness is simply given by that of the 
t-J model (up to a rescaling). Similarly the reduction of 
the spin stiffness due to the addition of impurities can be 
understood. 

Finally, to study the effects of the impurity on the ex- 
citations in the model we have computed the phase shifts 
acquired by holons and spinons due to scattering off the 
impurity in the case of a microscopic number of holes in 
the half filled ground state at vanishing magnetic field. 
The basic ingredient for this calculation is the quantiza- 
tion condition for factorized scattering of two particles 
with rapidities Ai and A2 on a ring of length N, namely 
exp(i7Vfc(Ai))S*(Ai - A 2 ) exp(-0(Ai)) = 1 where fc(A) is 
the physical momentum in the infinite periodic system, 
S(X) is the bulk scattering matrix for scattering of parti- 
cles 1 and 2 and ip(^i) is f ne phase shift acquired when 
scattering off the impurity (note that this incorporates 
the fact that there is no backscattering at the impurity) . 
Using the known result for S pl[ one extracts the impu- 
rity phase shifts using the method of [ p2|j2^ ]: both the 



spinon and holon impurity phase shifts are proportional 
to exp(— ik), where k denotes the physical momentum 
of spinons/holons in the t-J model without impurity. 
This reflects the fact that the impurity essentially de- 
couples from the chain at half filling leading to a chain of 
N — 1 sites. In addition the holons pick up a phase shift 
(2i\ — a)/(2iX + ia) due to the fact that the impurity 
site is charged. The pole at A — ia/2 corresponds to an 
impurity bound state for a < 2. 

To summarize, we have studied the effects of the addi- 
tion of integrable impurities to the supersymmetric t-J 
model on certain zero and finite temperature properties 
of the system. The properties of the impurity, which cou- 
ples to both spin and charge degrees of freedom, can be 
tuned by adjusting a continous parameter a. Compared 
to the "weak-link" type impurities investigated by Kane 
and Fisher [0 it appears to be very special in that its 
dc-conductivity is unchanged by the addition of a single 
impurity. We have argued that this is due to the absence 
of backscattering terms on the level of the dressed exci- 
tations (holons and spinons). Although the verification 
by explicit construction of the continuum limit for this 
system appears difficult one may speculate that similar 
to the case of a Kondo impurity in a Luttinger liquid || 
a backscattering term would drive the system to a new 
fixed point. Hence the present model can be interpreted 
as an unstable fixed point from a renormalization group 
point of view. 
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FIG. 2. (a) Impurity magnetization as a function of mag- 
netic field for band filling n e = 0.25 and several values of a. 
(b) Number of electrons located at the impurity as a function 
of the bulk electron-density for fixed magnetic field H — 0.1 
and several values of a. The dotted line denotes the critical 
electron-density n c below which the ground state is ferromag- 
netic. 




FIG. 1. Coupling of the impurity site (square) to the t-J 
bulk sites (circles). 
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FIG. 3. Charge-stiffness for a chain with 20 percent im- 
purities as a function of the electron density at H — 0.1 for 
several values of a. The dashed line denotes the stiffness at 
the critical electron density n c (a). (For comparison we have 
included the stiffness for the pure t-J chain D\ P j .) 



4 



